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Abstract 

We have studied whether the rotation and magnetic fields in neutrino-driven winds can be key processes 

for the rapid-process (r-process) nucleosynthesis. We have examined the features of a steady and subsonic 

wind solutions which extend the model of Weber and Davis (1967), which is a representative solar wind 

model. As a result, we found that the entropy per baryon becomes lower and the dynamical timescale 

becomes longer as the angular velocity becomes higher. These results are inappropriate for the production 

of the r-process nuclei. As for the effects of magnetic fields, we found that a solution as a steady wind 

from the surface of the proto-neutron star can not be obtained when the strength of the magnetic field 

becomes > 10 11 G. Since the magnetic field in normal pulsars is of order 10 12 G, a steady wind solution 

might not be realized there, which means that the models in this study may not be adopted for normal 

proto-neutron stars. In this situation, we have little choice but to conclude that it is difficult to realize a 

successful r-process nucleosynthesis in the wind models in this framework. 
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1. Introduction 
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It is one of the most important astrophysical problems that the sites where the rapid-process (r-process) nucle- 
osynthesis occurs are not still known exactly. There are, at least, three reasons that make the study on r-process 
nucleosynthesis important. One of them is a very pure scientific interest. The mass numbers of the products of 
r-process nucleosynthesis are very high (A = 80-250), which means that the most massive nuclei in the universe 
are synthesized through the r-process. You can guess easily that the situation in which the r-process nuclei are 
synthesized is a very peculiar one in the universe. We want to know where, when, and how the r-process nuclei are 
formed. Second reason is that some r-process nuclei can be used as chronometers. For example, the half-lives of 
232 Th and 238 U are 1.405 xlO 10 yr and 4.468 xlO 9 yr, respectively. Therefore, if we can predict the mass-spectrum 
of the products of r-process nucleosynthesis precisely, we can estimate the ages of metal-poor objects which contain 
the r-process nuclei by observing its abundance ratio. Third reason is that some r-process nuclei can be used as 
tools of the study on the chemical evolution in our Galaxy (e.g., Ishimaru and Wanajo 1999), which has a potential 
to reveal the history of the evolution of our Galaxy itself. Due to the reasons mentioned above, the study on the 
r-process nucleosynthesis is very important. 

The conditions in which the r-process nucleosynthesis occurs successfully are (e.g., Hoffman ct al. 1997): (i) 
neutron-rich (n n > 10 20 cm -3 ), (ii) high entropy per baryon, (iii) small dynamical timescale, and (iv) small Y e . 
This is because r-process nuclei are synthesized through the non-equilibrium process of the rapid neutron capture 
on the seed nuclei that is synthesized through the alpha-rich freezeout (e.g., Hoffman et al. 1997). In other words, 
an explosive and neutron-rich site with high entropy will be a candidate for the location where the r-process 
nucleosynthesis occurs. 

The candidates of the reliable sites are collapse-driven supernovae (e.g., Woosley et al. 1994) and/or neutron star 
mergers (e.g., Freiburghaus et al. 1999). This is because these candidates are thought to have a potential to satisfy 
the conditions mentioned above. However, we think that the collapse-driven supernovae are thought to be more 
probable sites than the neutron star mergers, because metal poor stars already contain the r-process nuclei (e.g., 
Freiburghaus et al. 1999). In fact, Mc William et al. (1995) reported that the abundance of Eu can be estimated in 
11 stars out of 33 metal-poor stars. These observations prove that r-process nuclei are produced from the early stage 
of the star formation in our Galaxy. Comparing the event rate of collapse-driven supernovae (10 -2 yr -1 Gal -1 ; 
van den Bcrgh and Tammann 1991) with the neutron star merger (10~ 5 yr -1 Gal -1 ; van den Heuvel and Lorimer 
1996; Bethe and Brown 1998), we can see that collapse-driven supernovae are favored since they can supply the 
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r-process nuclei from the early stage of the star formation in our Galaxy. Also, Cowan et al (1999) reported that 

the abundance ratio of r-process nuclei in metal poor stars are very similar to that in the solar system. This proves 

that r-process nuclei are synthesized through the similar conditions. This will be translated that, at least, most of 

the r-process nuclei are from one candidate. Therefore, we assume in this paper that most of the r-process nuclei 

are synthesized in the collapse-driven supernovae. 

There are many excellent and precise analytic and/or numerical computations on the r-process nucleosynthesis 
in the collapse-driven supernovae. However, so far it seems that there is no report that the r-process nuclei can be 
reproduced completely. For example, Takahashi et al. (1994) performed numerical simulations assuming Newtonian 
gravity and reported that entropy per baryon in the hot bubble is about 5 times smaller than the required value. 
Qian and Woosley (1996; hereafter QW96) also reported analytic treatments of the neutrino-driven winds from 
the surface of the proto-neutron star. At the same time, their analytical treatments are tested and confirmed by 
numerical methods. However, it was shown that the entropy derived by their wind solutions fall short, by a factor 
of 2-3, of the value required to produce a strong r-process (Hoffman et al. 1997). In order to solve this difficulty, 
Qian and Woosley (1996) included a first post-Newtonian correction to the equation of the gravitational force. As a 
result, they reported that the entropy increases and the dynamical timescale is reduced by a factor of ~ 2. Cardall 
and Fuller (1997) developed this argument by considering a fully general relativistic treatment. They showed that a 
more compact neutron star leads to higher entropy and a shorter dynamical timescale in the neutrino-driven wind. 
In order to confirm their conclusion quantitatively, Otsuki et al. (2000) have surveyed the effects of general relativity 
parametrically. They reported that r-process can occur in the strong neutrino-driven winds (L v ~ 10^ erg s ^) as 
long as a massive (^2.0 Mq) and compact (~ 10 km) proto-neutron star is formed. It is very interesting because 
such a solution can not be found in the frame work of Newtonian gravity (Qian and Woosley 1996). Such a solution 
is confirmed by the excellent numerical calculations (Sumiyoshi ct al. 1999). However, the equation of state (EOS) 
of the nuclear matter has to be very soft to achieve such conditions. Although a few non-standard models of EOS 
can satisfy them (Wiringa et al. 1988) as long as the matter is sufficiently cold, it seems to be very difficult to achieve 
them in the phase of the proto-neutron star. In fact, the r-process nuclei can not be produced in the numerical 
simulations with a normal EOS (Sumiyoshi et al. 1999). Thus, it seems that the difficulty can not be solved by only 
the effects of general relativity. 

There is only one report that r-process nucleosynthesis occurred successfully. That is the work done by Woosley 
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et al (1994; here after WWMHM94). In their numerical simulation, the entropy per baryon becomes higher and 

higher as the computation time goes on. Finally, at very late phase of neutrino-driven wind (~ 10 s after the 

core-collapse), successful r-process occurs. However, there are some problems in their results. First of all, it is 

unclear why the entropy per baryon at the late phase becomes so high as their results. In fact, when we adopt the 

analytic formulation of QW96, such a high entropy should not be obtained. Although the general relativistic effects 

are included in WWMHM94, such a high entropy could not be obtained in Otsuki et al (2000). Therefore, the 

discrepancy between WWMHM94 and QW96 can not be simply explained by only the general relativistic effects. 

Also, WWMHM94 has a problem that much nuclei whose mass numbers are ~ 90 are produced in the early stage 

of the neutrino-driven winds. To agree with the observational solar system abundances, we have to abandon the 

products at the early stage of the neutrino-driven winds. In addition, the successful mass-spectrum at the late 

phase of neutrino-driven winds would be destroyed when the reactions of neutral-current neutrino spallations of 

nucleons from He are taken into consideration (Meyer 1995). They reported that the entropy should be increased 

by (30-50)% in order to restore the A = 195 peak. They are extremely large modifications to the model. Although 

WWMHM94 is surely the very remarkable and interesting work, the problem of r-process nucleosynthesis has not 

been solved completely. 

Due to the reason mentioned above, it will be natural to think that there may be an (some) effect(s) that will 
help the r-process nucleosynthesis. In this study, we investigate the effects of rotation and magnetic fields on the 
synthesis of the r-process nuclei. In general, it is difficult to study their effects since the system including them is 
complicated and numerical simulations are needed in order to investigate them precisely. Since there is no numerical 
simulation like that, our final goal is to perform such realistic numerical simulations. However, as stated above, to 
perform such a numerical simulation will be a heavy task. Even if we can do it in future, it will be difficult to explain 
the results without any simple analytical studies. In this situation, before performing such numerical simulations, 
we examine the physical conditions of simple, exact, and steady solutions of the neutrino-driven wind including the 
effects of rotation and magnetic fields. We use the model that is the extension of the solution presented by Weber 
and Davis (1967), which is used as a representative model for the solar wind. In this study, we add the effects of 
neutrino heating and cooling to the solution and examine whether the effects of rotation and magnetic fields can 
help the synthesis of r-process nuclei. 

In section 5], we explain the formulation for the wind in the hot bubble. Results are shown in section EH. Summary 



No. ] Rapid-Process Nucleosynthesis in Neutrino-Magneto-Centrifugally Driven Winds 

and discussions are presented in section EJ. 
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Formulations 



2.1. Basic Equations 



In Gaussian units, the Euler equation acted on by electromagnetic forces can be written as (Shapiro and Teukolsky 
1983) 



^ - —VP - V$ - — VP 2 + — (B ■ V)P. 

at p Sup Anp 



Here 
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v • V 
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(2) 



alt dt 

is the Lagrangian time derivative following a fluid clement. 

In this paper, we study a steady flow which has ^-symmetry around the equatorial plane of the proto-neutron 
star. Thus, we use the spherical coordinate (r, 8, <f>) for convenience. The origin r=0 is set at the center of the 
proto-neutron star. In this coordinate, the Euler equation in the radial direction for the system that has (^-symmetry 
can be written as 



dv r 
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Here we used the conservation of mass: 
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The equation for the evolution of material energy, e, is 
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where q is the net specific heating rate due to neutrino interactions (Qian and Woosley 1996). In this study, we 
consider three neutrino heating and/or cooling processes (neutrino absorption on free nucleons, neutrino scattering 
processes on the electrons and positrons, and electron and positron capture on free nucleons) as 



<1 = QvN + <lve — 9eJV, 



(7) 
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where 

1 — x 

q vN = 1.55 x 10~ 5 N A [(1 - Y e )L VeM el eMeV + Y e L^ 51 el eMeV ] [erg s" 1 g" 1 ] , (8) 

T 4 / 6 \ 1 — x 

q ue = 3.48 x 10~ 6 7Va^^v L^ e;5iei/e; MeV + ip e .5i£p e .McV + =i^,5i£^,MeV —&r [ cr S ^ S^ 1 ] i ( 9 ) 

P8 V 7 / H 6 

and 

<7 e7V = 3.63 x IO-^aIScv [erg s" 1 g" 1 ] . (10) 

Here is the neutrino sphere radius in units of 10 6 cm, p s is the density in units of 10 8 g cm -3 , Tmcv is the 
temperature in units of 1 MeV, x = (1 — R 2 /z 2 ) 1 / 2 , Na is Avogadro number, L v ^\ is the individual neutrino 
luminosity in 10 51 erg s _1 , e^McV is an appropriate neutrino energy e„ in MeV (Qian and Woosley 1996). In 
this study, we set 4=0 at T <0.5 MeV, because free nucleons are bound into a-particles and heavier nuclei and 
electron-positron pairs annihilate into photons. 

The pressure P and internal energy e are determined approximately by the relativistic electrons and positrons 
and photon radiation as long as T > 0.5 MeV. Then, the pressure and internal energy can be written as 

IItt 2 fc 4 T 4 
180 h 3 c 3 

and 

117T 2 fc 4 T 4 



[dyncm^ 2 ] (11) 



f = 



[erg g- 1 ] , (12) 



60 n 3 ^p 

where k and h are Boltzmann and Planck constants, respectively. These are the basic equations in this paper. 
Precisely, although we might have to consider the effects of annihilation of electron-positron pairs on the dynamics 
at T < 0.5 MeV, the effects seems to be little (Sumiyoshi ct al. 1999). 

2.2. Model for the Wind 



In this study, we use the wind model presented by Weber and Davis (1967) which is the steady flow in the 
equatorial plane and has a velocity 

v = (v r ,0,v^) (13) 
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and a magnetic field 

B = (B r ,0,B 4 ). (14) 
We also assume that the system is symmetric with respect to the equatorial plane. In this case, Eq. (||) becomes 



dv r IdP GM vl Rx d 



•iv 



dr p dr r 2 r Aitpr dr 
Conservation of mass becomes 



(r£g. (15) 



de _ P_dp 
dr p 2 dr 



pv r r = f = const. (16) 
The equation of the evolution of material energy becomes 

(17) 

Furthermore, since div B = 0, 

r 2 B r = $ = const. (18) 
We also obtain from ^roti?^ = and Maxwell equation that 

r(t; r 5 - v^Br) = -VLr 2 B ri (19) 
where fl is the angular velocity of the surface of the proto-neutron star. From the steady-state 0-equation of motion, 

Br 

rVfj, ri?^ = L = const. (20) 

4tt pijy 

These equations in this subsection represent the steady neutrino-driven wind with proper magnetic fields in the 
equatorial plane that has ^-symmetry. See detail Weber and Davis (1967). 

2.3. Boundary Conditions 



In this study, the surface of the proto-neutron star is considered as the inner boundary. The inner boundary 
conditions are composed of density, luminosities of neutrinos, mass and radius of the proto-neutron star, velocity of 
the outflow, angular velocity, and strength of B r at the surface of the proto-neutron star. Temperature and electron 
fraction at the time when alpha-rich freezeout takes place are determined by these parameters as (Qian and Woosley 
1996) 

4,61^M.,MeV M (21) 
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(22) 



where L V fi\ is the individual neutrino luminosity in 10 51 ergs s" 1 , A = 1.293 MeV is the neutron-proton mass 
difference, and e^MeV is a neutrino energy in MeV. We assume that the neutron star radius is equal to the neutrino 
sphere radius. 

In this study, we assume that the luminosities of neutrinos are same (Qian and Woosley 1996; Otsuki et al. 2000). 
The energy of neutrinos are assumed to be 12, 22, and 34 MeV for v e , z? e , and other neutrinos, respectively (Woosley 
et al. 1994; Qian and Woosley 1996; Otsuki et al. 2000). Surface density is assumed to be 10 10 g cm -3 (Otsuki et 
al. 2000). In the previous works, initial velocity of the outflow is chosen so that tvo, which is the radial velocity at 
the surface of the proto-neutron star, becomes less than t> r o,crit- ^ro.crit is the critical value for supersonic solution 
(Qian and Woosley 1996; Otsuki et al. 2000). In this work, we also adopt this assumption so that the flow becomes 
subsonic and contains no critical point. In particular, we take v t q as ivo ~ i>ro,crit i n this study. This means that the 
initial velocity (velocity at the surface of the proto-neutron star) is set to be maximal one because surface density 
is set to be constant (10 10 g cm -3 ). In case wc try to survey the flow that contains a transition point like a shock 
front, we can not use Eqs. (|l^)~ (|2^)- This is because these differential equations diverge and break down. In order 
to treat such a flow that contains a discontinuity, we have to use the Rankine-Hugoniot relation instead of these 
equations. We will examine such flows in the forthcoming paper. Other parameters are changed parametrically. 
The parameters and the name of the models are given in Table 1. We take r — 10 9 cm for the radius of the outer 
boundary (Otsuki et al. 2000). Since the radius of the Fe core is about 10 s cm, we think the radius of the outer 
boundary is large enough to investigate the r-proccss nucleosynthesis that occurs in the hot bubble. 



Output parameters are shown in Table 2. Entropy per baryon (£), dynamical timescale (rdyn), electron fraction 
(Y e ), and temperature (T) at the outer boundary (r = 10 9 cm) are shown in the table. As long as T > 0.5 MeV, 
entropy per baryon in radiation dominated gas is written by 



3. 



Results 
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(23) 
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where tun is the nucleon rest mass. Even if electron-positron pairs are disappeared for T < 0.5 MeV, the entropy per 

baryon is conserved because the wind expands adiabatically in this region. We define the dynamical timescale to be 

the elapsed time for the temperature to decrease from 0.5 MeV to 0.2 MeV. This is because r-process nucleosynthesis 

occurs in this temperature range (Woosley et al. 1994; Takahashi et al. 1994; Qian and Woosley 1996). In the 

following subsections, the influences of the luminosity of neutrinos, rotation, and magnetic fields on the dynamics 

are discussed, respectively. 

3.1. Influence of luminosity of neutrinos 



In order to examine the influence of the luminosities of neutrinos on the dynamics, we focus on the results of 
Models lOAa, lOBa, and lOCa. When the angular velocity and magnetic fields are sufficiently small, the solution 
have to agree with QW96. Then, the results of Models lOAa, lOBa, and lOCa should be explained by the solution 
obtained in QW96. In QW96, the entropy per baryon at the beginning of the alpha-rich freezeout, dynamical 
timescale, and radial velocity at the surface of the proto-neutron star are estimated as 

S/k ~ 235L^ 51 e p e W#6 3 ( TTTT~ ) > ( 24 ) 



1.4M 



and 



4 /I0 10 gcur 3 \ /1.4M(T X 2 



v rQ ~ 1.8Ll e)51 ep 3 e>M eV^ * [ ~ p ) [-Jf^ ) K s" 1 ] , (26) 

respectively. 



From Eqs. (24)-([26j) and Table 2, we can see that the dependence of the entropy per baryon, dynamical timescale, 
and initial velocity on the luminosity of neutrinos arc well explained by these equations. We also note that these 
values in Table 2 are well reproduced by these equations within a factor of 2-3, which warrants the accuracy of our 
calculations. 



3. 2. Influence of rotation 
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The effects of rotation in this framework can be understood by comparing the results of Models lOAa-lOAc, 

lOBa-lOBd, and lOCa-lOCd in Table 2. As is clear from the results of Models lOAc, lOBd, and lOCd, in which the 

rotation period of the proto-neutron star is (0.5 - 1) ms, the entropy per baryon becomes lower and the dynamical 

timescale becomes longer as the angular velocity becomes higher. This tendency is inappropriate for the success of 

the r-process nucleosynthesis (see section 1). 

We consider the reason for this tendency. At first, we show in Figure 1 the outflow velocity, temperature, and 

density as a function of r for Model lOBa and Model lOBd. Entropy per baryon as a function of radius (r) from the 

center of the proto-neutron star in these models is also shown in Figure 2. From these figures, we can see that the 

density in Model lOBd is much higher than that in Model lOBa at relatively small radius (r < 5 x 10 7 cm), which 

results in the lower entropy per baryon (see also Eq. (p3|)). From Eq. (|l6|), we can see that v r does not increase so 

much if the density does not decrease, which is verified in Figure 2. Since Td yn becomes longer when v r is slower 

in the range T = (0.2 - 0.5)MeV, the dynamical timescale in Model lOBd is longer than that in Model lOBa. This 

is the reason for the tendency mentioned above. All we have to do is to find the reason why the density in Model 

lOBd is higher than that in Model lOBa at small radius. 

From the basic equations in subsection 2.2, the derivatives of p by r can be written as 

if _ p± _ gm p , 2/$n P a^, pv% 

dp pr 5 ev T r 2 47r/ 2 — p$ 2 r 



(27) 

dr A^P^L-nr^B^ P_ ( p\_ p v ' 

r(47r/ 2 — p$ 2 ) 2 ' tp \ ' ' pi p 2 r l 

We show the absolute value (in cgs units) of the each component in Eq. ( p7| ) for Model lOBa in Figure 3a. From 
this figure, we find out what component dominantly contributes to the gradients of density. To compare them, 
we also show those for Model lOBd in Figure 3b. Lines (a)-(f) correspond to 2/ 2 /pr 5 , Pq/ev r , GMp/r 2 , pvi/r, 
P(e + P/p)/ep, and f 2 p 2 r A as a function of r, respectively. Here we note that is set to be zero in these models. 
Therefore, there is no components that include B^ in the figures. Line (d), which represents pv^/r, is very important 
here. As can be seen from Figure 3, the density gradient in both models can be approximated by 

, _GMp P4 

J r> + I . (28) 

dr 

Since the value for pv^/r is larger in Model lOBd than in Model lOBa, the absolute value of the density gradient in 
Model lOBd is smaller than in Model lOBa. Then, the density in Model lOBd becomes larger than that in Model 
lOBa for a relatively small r because the density at the surface of the proto-neutron star is set to be 10 10 g cm~ 3 in 
both models. This is the reason for the tendency mentioned above. 
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3. 3. Influence of luminosity of magnetic field 
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The effects of magnetic fields in this framework can be understood by considering the results of Models lOAe, 
lOBg, and lOCf in Table 2. In Models lOAe, lOBg, and lOCf, in which the strength of the magnetic field is ~ 10 11 
G, we could not find a solution as a steady wind from the surface of the proto-neutron star because of the following 
reasons. First, we focus on the first term of the denominator in Eq. (|27|). It diverges when 47r/ 2 — p$ 2 becomes 
zero. When this relation is satisfied, the solution diverges. This condition can be rewritten as 

' 10 10 g cm" 3 \ 2 / 10 5 cm s" 1 \ 2 / B r0 



1 - 7.96 x 10 



-8 



0, (29) 



Po J V v r0 ) V10 12 G, 
where poj v r o, and B r o are the density, radial velocity, and radial magnetic field at the surface of the proto-neutron 
star. 

From Eq. (^9|), we can see that the solution does not diverge when B r0 is small enough (Models in this study 
except for Models lOAe, lOBg, and lOCf). This is because the second term in Eq. ( p9|) is small enough. In Models 
lOAe, lOBg, and lOCf, initial velocities have to be set to be high in order to avoid the divergence. However, when 
the initial velocity is set to be so high that the second term in Eq. ( p9[) becomes negligible, the solution can not 
satisfy the condition that the solution is subsonic. That is, v r o becomes larger than u r o,crit an d the value of the 
denominator of Eq. ( p7| ) becomes zero at the critical point and the solution diverges after all (see also QW96). This 
is the reason why a steady solution can not be obtained when the magnetic fields becomes strong enough. It is also 
noted that we can not set tvo to be small enough to avoid the divergence resulting from Eq. ( p9[ ) , because the initial 
velocity is approximated well by Eq. ([26]), which is confirmed by precise numerical simulations (Qian and Woosley 
1996; Sumiyoshi et al. 1999). 

4. Summary and Discussion 



We have studied whether the effects of rotation and magnetic fields on the neutrino-driven winds could be necessary 
ones for the r-process nucleosynthesis. We have studied the effects of the rotation and magnetic field using simple 
models which are the extensions of the solution presented by Weber and Davis (1967), because the results of a realistic 
numerical simulation concerning with such a topic will not be understood clearly without any analytical studies. 
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Although our final goal is to perform such realistic numerical simulations, this approach would be a necessary step 

to understand the effects of the rotation and magnetic fields on the r-process nucleosynthesis. 

As a result, we found that the entropy per baryon becomes lower and the dynamical timescale becomes longer as 

the angular velocity becomes higher, which is a bad tendency for the success of the r-process nucleosynthesis. This 

is because the absolute value of the density gradient becomes smaller due to the effects of rotation and the density 

is kept to be high at relatively small radius (r < 5 x 10 7 cm), which results in the lower entropy per baryon and 

longer dynamical timescale. As for the effects of magnetic fields, we found that a solution as a steady wind from 

the surface of the proto-neutron star can not be obtained when the strength of the magnetic field becomes > 10 11 

G. This is because the density gradient diverges at the critical point, which emerges at low density region like the 

circumstances of neutrino-driven winds in this study as long as the amplitude of the magnetic field is large enough. 

As a conclusion, we have to say that it seems to be difficult to cause a successful r-process nucleosynthesis in the 

wind models in this study. 

Since the magnetic field in normal pulsars is of order 10 12 G (e.g., Thompson 2001), the fact that a steady wind 
solution can not be obtained as long as the radial component of the magnetic fields at the surface of the proto-neutron 
star is larger than 10 11 G seems to mean that the models in this study may not be able to be adopted in many cases. 
However, our models could be used since radius of proto-neutron stars are tend to be larger than normal pulsars 
(e.g., Wilson 1985) and the resulting magnetic fields are weaker. 

We have to emphasize that there are some assumptions in this study. So, we can not say that we have proved 
that a successful r-process nucleosynthesis does not occur in neutrino-driven winds in which the effects of rotation 
and magnetic fields are taken into consideration. For example, we assumed that the flow is subsonic and there is no 
critical point, which is the common assumption in the previous studies on the r-process nucleosynthesis (Qian and 
Woosley 1996; Otsuki et al. 2000). However, we think that we do not need to restrict the solutions in such a way, 
that is, there may be a transition point at which Eqs. (p^)-(|20|) break down. In most cases, the transition point 
will be a shock front. It means that the flow will gain entropy at the transition point, which will be a good sense to 
produce the r-process nuclei. The problem whether the flow contains transition points or not depends sensitively on 
the initial velocity on the surface of the proto-neutron star. So our final goal is to determine physically the velocity 
at the surface of the proto-neutron star. It means that the tvo should not be given as an input parameter. It should 
be an output parameter. We have to investigate the mechanism for determining the outflow velocity at the surface 
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of the proto-neutron star for further discussions. We also assumed that the flow is steady. We should also investigate 

the features of the unsteady flows, although the flow is assumed to be steady in this study. It will be investigated 

by numerical tests assuming a simple environment. We will perform such numerical tests in the near future. 

We have assumed that the form of the wind is similar as the solution of Weber and Davis (1967) here. Of course, 

there will be a variety of flows in which the effects of rotation and magnetic fields arc taken into consideration, 

including the jets (Nagataki 2001). So it will be worth while surveying physical conditions using a variety of forms 

of the neutrino-driven winds, because we will be able to understand more precisely and correctly the results of the 

realistic numerical simulations. Of course, our final goal is to perform realistic numerical simulations for the wind 

in a collapse-driven supernova and seek a final answer whether the rotation and magnetic fields in neutrino-driven 

winds are the key processes for the r-process nucleosynthesis or not. 

This research has been supported in part by a Grant-in-Aid for the Center-of-Excellence (COE) Research 
(07CE2002) and for the Scientific Research Fund (7449, 199908802) of the Ministry of Education, Science, Sports 
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Table 1. Model Names and Input Parameters. 





Mass 


Radius 




v r0 


n/2ir 


B r0 


Model 


(M Q ) 


(km) 


(10 51 ergs s- 1 ) 


(cm s _1 ) 


(Hz) 


(gauss) 


lOAa 


1.4 


10 


3.0 


8.30(+5) 


1.0(+1) 





lOAb 


1.4 


10 


3.0 


8.30(+5) 


1.0(+2) 





10 Ac 


1.4 


10 


3.0 


1.10(+6) 


1.0(+3) 





10 Ad 


1.4 


10 


3.0 


8.30(+5) 


1.0(+1) 


5.0(+10) 


lOAe 


1.4 


10 


3.0 


8.50(+5) 


1.0(+1) 


5.0(+ll) 


lOBa 


1.4 


10 


1.0 


1.24(+5) 


1.0(+1) 





lOBb 


1.4 


10 


1.0 


1.24(+5) 


1.0(+2) 





lOBc 


1.4 


10 


1.0 


1.66(+5) 


1.0(+3) 





lOBd 


1.4 


10 


1.0 


6.30(+5) 


2.0(+3) 





lOBc 


1.4 


10 


1.0 


1.24(+5) 


1.0(+1) 


1.1(+11) 


lOBf 


1.4 


10 


1.0 


1.25(+5) 


1.0(+1) 


1.2(+11) 


lOBg 


1.4 


10 


1.0 


1.20(+5) 


1.0(+1) 


1.3(+11) 


lOCa 


1.4 


10 


0.6 


5.20(+4) 


1.0(+1) 





lOCb 


1.4 


10 


0.6 


5.20(+4) 


1.0(+2) 





lOCc 


1.4 


10 


0.6 


7.00(+4) 


1.0(+3) 





lOCd 


1.4 


10 


0.6 


2.70(+5) 


2.0(+3) 





lOCe 


1.4 


10 


0.6 


5.20(+4) 


1.0(+1) 


5.0(+10) 


lOCf 


1.4 


10 


0.6 


1.00(+4) 


1.0(+1) 


6.0(+10) 



Model names and input parameters. Mass and radius of the proto-neutron star, luminosity of the anti-electron 
neutrino, radial and angular velocity, and radial component of the magnetic field at the surface of the proto-neutron 
star, respectively. 
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Table 2. Model Names and Output Parameters. 





c 
o 


T dyn 




+ b 


A !X s~\ r\ i~\ 1 

iviociei 


(k) 






Me V 


1 n A n 
lUAa 


oz 


o o/ oa 

Z.Z(-Z ) 


U.4o 


a a ( oA 
0.4(-z J 


in AV. 
1UAD 


QQ 

oo 


O Q / OA 

z.o(-z j 


U.4o 


7 ft / oA 
|\0(-Z J 


1UAC 


to 


r r\( oA 

o.U(-z J 


U.4o 


no/ OA 

y.o(-z J 


lUAa 


QO 
oz 


O ( OA 

z.z^-z ) 


U.4o 


ft /I / OA 
0.4(-z J 


1 n A 
lUAe 






U.4o 




ll)r>a 


no 

yy 


K K 1 OA 

o.o(-z ; 


U.4o 


A 7/ OA 

4. / (-z J 


1U-DD 


yy 


O.O^-z ) 


O /I Q 
U.4o 


/I /I f OA 
4.4^-z j 


lUrSC 


QQ 
OO 


K K ( OA 

o.o^-z J 


U.4o 


Q ft/ OA 
o.O^-z j 


1U-DC1 


OO 


7 7/ OA 
' • ' 1"^ ) 


U.4o 


O O / OA 

Z.Z\-Z ) 


lUrSC 


on 

yy 


K Kf OA 
O.O^-Z J 


U.4o 


A 7/ OA 

4. I (^-z j 


1U-DI 


y ( 


o n/ oA 

o.y (-z j 


U.4o 


o n./ oA 
z.LH-z J 


lOtfg 






U.4o 




lOCa 


107 


8.7(-2) 


0.43 


4.1(-2) 


lOCb 


108 


9.9(-2) 


0.43 


5.1(-2) 


lOCc 


94 


6.9(-2) 


0.43 


1.2(-2) 


lOCd 


58 


1.2(-1) 


0.43 


1.5(-2) 


lOCe 


107 


8.7(-2) 


0.43 


4.1 (-2) 


lOCf 






0.43 





Name and output parameters of Models. Entropy per baryon, dynamical timescale (rd yn ), electron fraction, tem- 
perature at the outer boundary are shown respectively. The reason why physical quanta are not written in some 
models is that steady solutions can not be obtained in these models. 



Figure Captions 

Fig. 1. Outflow velocity, temperature, and density as a function of radius (r) from the center of the proto- neutron 
star. These values are written in unit of 10 7 cm s _1 , 1 MeV, and 10 8 g cm -3 , respectively. Solid lines correspond 
to Model lOBa, whereas dashed lines correspond to Model lOBd. 

Fig. 2. Entropy per baryon as a function of radius (r) from the center of the proto-neutron star. Solid line 
corresponds to that of Model lOBa. Dashed line corresponds to that of Model lOBd. 

Fig. 3. Upper panel: absolute value (in cgs units) of the each component of Eq. (27) for Model lOBa. Lines (a)-(f) 
correspond to 2/ 2 /pr 5 , Pq/ev r , GMp/r 2 , pv^/r, P(e + P/p)/cp, and f 2 p 2 r 4 as a function of r, respectively. 
The discontinuity of line (b) at r ~ 10 7 cm reflects the frcczcout of the neutrino reactions (see subsection 2.1). 
Lower panel: same as left, but for Model lOBd. 
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